Astronomy Letters, Vol. 33, No. 5, 2007, pp. 299-308. Translated from Pis'ma v Astronomicheskii 
Zhurnal, Vol. 33, No. 5, 2007, pp. 340-351. 



High-Mass X-ray Binaries and the Spiral Structure of the Host 

Galaxy 



P. E. Shtykovskiyi 2*, M. R. Gilfanov^-^ 

^ Space Research Institute, Profsoyuznaya str. 84/32, Moscow 117997, Russia 
^MPI for Astrophysik, Karl-Schwarzs child str. 1, Garching, 85741, Germany 



We investigate the manifestation of the spiral structure in the distribution of high-mass 
X-ray binaries (HMXBs) over the host galaxy. We construct the simplest kinematic model. 
It shows that the HMXBs should be displaced relative to the spiral structure observed in 
such traditional star formation rate indicators as the Hq, and far-infrared emissions because of 
their finite lifetimes. Using Chandra observations of M51, we have studied the distribution 
of X-ray sources relative to the spiral arms of this galaxy observed in Hq,. Based on K- 
band data and background source number counts, we have separated the contributions from 
high-mass and low-mass X-ray binaries and active galactic nuclei. In agreement with model 
predictions, the distribution of HMXBs is wider than that of bright HII regions concentrated 
in the region of ongoing star formation. However, the statistical significance of this result 
is low, as is the significance of the concentration of the total population of X-ray sources 
to the spiral arms. We also predict the distribution of HMXBs in our Galaxy in Galactic 
longitude. The distribution depends on the mean HMXB age and can differ significantly 
from the distributions of such young objects as ultracompact HII regions. 
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INTRODUCTION 



Chandra and XMM-Newton X-ray observations of nearby galaxies have revealed rich 
populations of compact X-ray sources in them (see, e.g., Fabbiano 2006; Kilgard et al. 2005). 
Their spectral analysis, flux variability, luminosity functions, and spatial distribution make 
it possible to identify them with known (in our Galaxy) classes of objects - high-mass and 
low-mass X-ray binaries, supernova remnants, ultrasoft X-ray sources, and to distinguish a 
hitherto unknown class of ultr aluminous X-ray sources. 

From the standpoint of their relation to star formation, the X-ray sources can be divided 
into two groups: young objects, such as high-mass X-ray binaries (HMXBs) and supernova 
remnants (SNRs), and old objects, such as low-mass X-ray binaries (LMXBs). It would be 
natural to expect the former to concentrate in regions of recent star formation and the lat- 
ter to follow the galactic stellar mass distribution. Indeed, Chandra observations of nearby 
galaxies show that the populations of X-ray sources in starforming galaxies differ radically 
from those in elliptical galaxies. This is particularly clearly seen when the luminosity func- 
tions of these sources are analyzed: the luminosity functions for the former are in the form 
of a power law with an index of 1.6, as for the HMXBs in our Galaxy, while those for the 
latter are similar in form to the luminosity function of LMXBs in our Galaxy (Grimm et al. 
2003; Gilfanov 2004). 

In nearby spiral galaxies, two components should be observed in the distribution of X-ray 
sources: HMXBs concentrating to the spiral arms and LMXBs distributed more smoothly 
with the maximum density at the galactic center. On the one hand, Chandra observations 
of such galaxies as M51 and MlOl suggest that the X-ray sources actually concentrate to 
the spiral arms. However, a more detailed analysis indicates that, in reality, the picture 
is more complex. First, HMXBs reflect the star formation in a galaxy that took place 
~ 5 — 60 Myr ago, i.e., strictly speaking, they are not an instantaneous star formation rate 
(SFR) indicator (Shtykovskiy and Gilfanov 2005). Such times may turn out to be important 
from the standpoint of galactic dynamics, since the characteristic revolution time of the stars 
in a typical spiral galaxy is ~ 200 Myr. As a result, the spiral structure observed in the 
HMXB distribution will be distorted and displaced relative to that observed, for example, 
in Hq. On the other hand, the density contrast between the arms and the interarm space 
in the stellar population of certain galaxies can reach ~ 3. This will give rise to a spiral 
structure in the old population of LMXBs as well. 

In this paper, we construct a kinematic model for the spatial distribution of X-ray binaries 
in a spiral galaxy, compare its predictions with Chandra observations of M51, and make 
predictions for our Galaxy. 

1. THE SPIRAL STRUCTURE IN VARIOUS SFR INDICATORS 

The basic principles of the modern theory of spiral structure in galaxies were laid down 
by Lin and Shu (1964, 1966) and Lin et al. (1969), who suggested the hypothesis of a 
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quasi-stationary density wave. According to this hypothesis, the spiral structure of a galaxy 
is a manifestation of the density wave - a stellar disk and gas density perturbation that 
propagates through the galaxy and that does not decay for a long period. Although the 
reasons why a stationary density wave emerges and is maintained are unknown, there are 
many candidates for their role, such as the influence of a neighboring galaxy, asymmetry at 
the galactic center, etc. Since the velocity dispersion in the interstellar gas is comparatively 
low, the amplitude of the density wave triggered by a gravitational potential perturbation 
can be large. In contrast, the amplitude of the density wave in the stellar disk should be 
small (at least for an isolated galaxy). Moreover, a shock wave leading to a very narrow 
zone of gas compression with a large density jump can generally emerge in the gas (Roberts 
1969). The gas compression triggers star formation in it; as a result, young stars producing 
a distinct galactic spiral structure in the optical band are formed. 

In addition to the visible range, the spiral structure can also be observed in many other 
SFR, stellar mass, and gas indicators. The H^, ultraviolet, and far-infrared emissions pro- 
duced by massive stars, the near-infrared emission from old stars, the 21-cm HI emission, and 
the CO emission are the most important ones. Simple considerations suggest that the spiral 
structure will be different in different indicators. For example, the 11^^ emission originates 
in HII regions containing young stars with masses 10 Mq, while the ultraviolet emission is 
associated with the photospheric emission of stars from a wide mass range (Kennicutt 1998). 
Since the angular velocities of the stellar disk and the density wave are different, young 
stars are displaced from the density wave shock front with time. Due to short lifetimes of 
the massive stars responsible for the emission, 20 Myr, its peak should lie not far from 
the spiral density wave front. On the other hand, the less massive stars responsible for the 
ultraviolet emission have lifetimes up to 100 Myr and can be displaced significantly from 
their birthplaces. As a result, the spiral structure in the ultraviolet will be considerably 
wider than that in H^. Displacements of this kind between various indicators were actually 
observed, for example, between Hq and the ultraviolet in M51 (Petit et al. 1996), HI and 
the nonthermal radio continuum inM51 and M83 (Tilanus and Allen 1991), etc. 

The situation with HMXBs is similar to that described above. Indeed, the lifetime of 
an optical counterpart with a mass of ~8 Mq is ~ 40 Myr, i.e., HMXBs reflect the star 
formation that took place several tens of Myr ago. A clear example of this, though unrelated 
to the spiral structure, is the absence of correlation between the surface density of HMXBs 
and the Ha emission intensity in the Large Magellanic Cloud (Shtykovskiy and Gilfanov 
2005). 

1.1. The Spiral Structure in the Distribution of High-Mass X-ray Biuciries 

To summarize the aforesaid, we conclude that the following factors affect the spatial 
distribution of HMXBs. 

(i) The spatial distribution of star-forming regions. The theory of density waves and ob- 
servations of spiral galaxies show that a logarithmic spiral can serve as a good approximation 
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for the location of the spiral density wave front. The location of a spiral arm is then specified 
in polar coordinates (r, 6) by the relation 



where ijj is the angle between the spiral and the tangent to the circumference called the pitch 
angle and (jq, ©q) is a point on the spiral, for example, its origin. We assume that the star 
formation takes place in a narrow region along this spiral. 

(ii) The dynamics of HMXBs. The HMXBs born at the shock front are drawn into the 
overall motion of the stellar disk in accordance with the galactic rotation curve. To a first 
approximation, it can be represented as the motion in circular orbits around the galactic 
center. The spiral density wave rotates with the pattern speed Qp in the same direction as 
the stellar disk. The HMXB locations relative to the instantaneous location of the spiral 
density wave front, which is initially equal to zero, in time r is then given by the following 
equation: 



where A© is the displacement of the objects relative to the instantaneous location of the 
density wave and Jl(r) is the galactic rotation curve. In other words, this formula gives 
the locations of objects with age r relative to the currently observed starforming region. 
Obviously, the HMXBs move faster than the density wave within the corotation radius 
{fl{rcr) — ^p) and lag behind it outside Tcr- 

In this formula, the gravitational field perturbation in the galaxy and the fact that the 
stars move in epicyclic orbits, not in circular ones, are ignored. To accurately describe the 

dynamics of the various components of a spiral galaxy, we must also consider the dynamics 
of gas clouds and the star formation in them under the effect of gravitational instability. 
Models of this interaction (see, e.g., Leisawitz and Bash 1982; Roberts and Stewart 1987) 
yield contradictory results: young stars will move almost along the spiral arm for a certain 
time in some of the models and will be rapidly displaced from it in other models. 

(iii) The evolution of the HMXB number with time elapsed since the star formation event. 
A detailed investigation of this question is a difficult task that requires developing population 
synthesis models and is beyond the scope of this paper. Here, we only briefly note the main 
factors that specify the evolution time scale of the population of HMXBs. First, this is the 
time it takes for a black hole or a neutron star to be formed. Here, two natural time scales 
can be identified: the hfetime of the most massive star (r^IOO Mq) corresponding to the 
formation of the first compact object (black hole), tmin ~ 2 — 3 Myr, and the lifetime of the 
least massive star («i8 Mq) capable of forming a compact object - the formation time of the 
last neutron star, t^nax ~ 40 Myr. The evolution of the companion star until the onset of an 
active phase, mass accretion onto the compact object, is of no less importance. The duration 
of this phase can also be rather long, up to several tens of Myr, depending on the mass of 
the companion star. In addition, the duration of the X-ray activity phase, which should 
be much shorter than the first two phases, 10^ — 10^ yr, must be taken into account. As a 
result, the maximum age for a HMXB can reach more than 40 Myr. This is also confirmed 



© — ©0 = ln(r/ro)/ tanip, 



(1) 



A© = (Q(r) - Qp)T, 



(2) 
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Fig. 1: Spiral structure in the distribution of HMXBs. Different curves corresponds 
to the locations of objects with ages of 10, 20, and 40 Myr calculated using Eq. (2). 
The location of the region of ongoing star formation (black curve) is also shown. The 
circumference indicates the location of the corotation radius. 

by observations: for example, analysis of the population of HMXBs in the Magellanic Clouds 
showed that binaries with ages up to '^60 Myr could be present among them (Shtykovskiy 
and Gilfanov 2007). Below, as the time scale for the existence of a HMXB population, we 
use both a conservative estimate of tmax ~ 40 Myr and larger values, demonstrating the 
magnitude of the effect in the latter case for the oldest binaries. 

It is also worth noting that the displacement of HMXBs relative to other SFR indicators 
observed in galaxies at distances ^ 1 Mpc can be somewhat smaller than that in nearby 
galaxies. Indeed, when X-ray binaries are observed in distant galaxies, there is a selection 
effect: the brightest objects are selected. Since the luminosity of a HMXB is related to 
the mass of its optical counterpart (Postnov 2003), a mass selection effect also arises, i.e., 
predominantly the highest-mass and, accordingly, youngest stars for which the displacement 
from the spiral arms is smaller, will be observed. However, apart from the mass of the 



- 5 - 




Fig. 2: (a) Distribution of distances from bright HII regions (Petit et al. 1996) to 
the nearest spiral arm (solid histogram) in M51. Also shown is the distribution for 
objects with an age of 40 Myr (dotted histogram), (b) The same for objects within 
the corotation radius. 



counterpart, the luminosity also depends on the orbital size. This reduces the selection 
effect and increases the contribution from "old" HMXBs to the observed population of X- 
ray sources. A detailed analysis of this effect is beyond the scope of this paper. 

To estimate the displacement of X-ray sources relative to the galactic spiral structure, we 
calculated the expected spatial distributions of HMXBs with ages of 10, 20, and 40 Myr. 
As the region of ongoing star formation, we used the approximation of the inner part of the 
star-forming region in M51 (see below) by a logarithmic spiral {ip = 19 °). The displacement 
of HMXBs relative to the instantaneous location of the spiral was calculated using Eq. (2) 
by assuming that the rotation curve and the density wave pattern speed also corresponded 
to those of M51. The derived distributions together with the current location of the spiral 
structure are shown in Fig. 1. 



2. THE SPIRAL STRUCTURE AND HIGH-MASS X-RAY BINARIES: 
COMPARISON WITH OBSERVATIONS 



Various methods can be used to study the spiral structure in the observed distribution 
of sources or emission intensity at a certain wavelength. For example, the distribution of 
displacements of the polar angle relative to the spiral angle, x = ~ ln(r/ro)/ tanip, can be 
constructed. In this presentation method, a logarithmic spiral with two arms is projected 
into two 5-functions near and vr. However, when the shape of the spiral structure is far 
from a logarithmic spiral, it is convenient to use the distributions of source distances to it. 
Obviously, this distribution will change in pattern as the sources "grow old" . This effect can 
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be demonstrated with an actual example using bright HII regions, for which the distribution 
of distances to the nearest spiral arm characterizes the locations of the youngest objects. The 
distribution for older objects can be obtained by displacing their locations relative to the 
current location of the spiral in accordance with Eq. (2). For this kind of demonstration, 
we chose bright HII regions in M51 ((S> 1.5 X 10"^^ W/m^) from the catalog by Petit 
et al. (1996) and displaced them by assuming that the spiral was a trailing one (i.e., the 
rotation is counterclockwise; see Fig. 3) . We then constructed the distributions of distances 
from the original and displaced sources to the nearest spiral arm (for more detail on the 
spiral structure in M51, see below) for sources in a wide range of galactocentric distances, 
2.3 kpc<r<9.5 kpc, and within the corotation radius, 2.3 kpc<r<5.6 kpc. A distortion of 
the distribution is clearly seen in both cases (see Fig. 2); in the latter case, it is in the form 
of a displacement. 

Obviously, for comparison with the qualitative model constructed above, a galaxy with 
active star formation and a low stellar mass whose spiral structure is close to a logarithmic 
spiral with a large pitch angle would be ideal. In addition, it is desirable that the angular 
size of the galaxy should be small enough for the contribution from background sources to 
be at a minimum. 

2.1. Compcirison with observations: M51 

One of the most suitable galaxies for our purpose among those observed by Chandra is 
M51. This galaxy has a distinct spiral structure and is located at a distance of 9.7 Mpc 
(Sandagc and Tammann 1975). The displacements of various indicators of the spiral struc- 
ture in M51 clearly indicate that we arc dealing with a density wave (sec. e.g., Tilanus and 
Allen 1989, 1991). To determine the orientation of M51 in space, we take the coordinates 
of its center RA(J2000)=13'*29"*52^71 and DEC(J2000)=47°25'42".6 (Ford et al. 1985), 
P.A. = 170°, and its inchnation to the plane of the sky i=20° (TuUy 1974). As the spiral 
wave pattern speed, we use Qp — 38 km/s/kpc obtained by Zimmer et al. (2004) by the 
Weinberg- Tremain method. Data on the rotation curve were taken from Tilanus and Allen 
(1991) and fitted by the law V{r) oc y^{r/roy--^/{l + (r/rg)^-^). The rotation curve and the 
density wave pattern speed correspond to a corotation radius icr ~ 5.6 kpc, in agreement 
with the values obtained by other methods (Vogel et al. 1993). 

An important peculiarity of M51 is its interaction with its companion, NGC 5195. As a 
result, it is often cited as an example of a galaxy where the spiral structure can be excited 
and maintained by tidal forces (see, e.g., Toomrc 1978). The tidal interaction is probably 
also responsible for the unusually high density contrast in the K band and, hence, in the 
stellar population of M51 (Rix and Rieke 1993). Obviously, a high density contrast in the 
stellar disk can change radically the manifestation of the spiral structure in X-ray binaries. 
For example, it would be natural to expect the appearance of a spiral structure in the 
distribution of not only HMXBs, but also LMXBs. 

The spiral structure in M51 is more complex in shape than a simple logarithmic spiral 
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and can be approximated by the latter only in a limited range of radii. Moreover, the ionized 
gas velocity field on the periphery of M51 indicates that the outermost parts of the spiral 
rotate with a pattern speed different from that of the inner spiral (Vogel et al. 1993). 

Below, we exclude the region around NGC 5195 with a radius of ^ 1.7' from our analysis. 

2.1.1. X-ray sources toward M51. 

For our analysis, we used the list of sources detected by Chandra (Terashima and Wilson 
2004). The total number of sources within a galactocentric radius of ^ 10 kpc is 88 at 
the limiting sensitivity S^m — 1-4 x lO"-*^^ erg/cm^/s in the 0.5-8 keV energy band (for a 
photon index of 1.5 and Nh — 1.3 x 10^° cm~^), corresponding to the luminosity Lx — 
1.6 X 10^^ erg/s. As with any other galaxy, most of the X-ray sources toward M51 belong to 
one of the following classes: HMXBs, LMXBs, SNRs, and background active galactic nuclei 
( AGNs) . Below, we discuss the contribution from each of them to the population of X-ray 
sources toward M51. 

(1) AGNs. According to the background source number counts by Moretti et al. (2003), 
^cxB ~13 cosmic X-ray background (CXB) sources are expected in a circle with a radius 
of 3.4' if we use the source counts in the soft 0.5-2 keV energy channel (S;^^ = 4.62 x 
20-16 erg/cm^/s, assuming a = 1.5 and Nh = 1-3 x 10^°cm~^) and ^cxb ~22 background 
sources are expected if we use the source counts in the hard 2-10 keV energy channel (S^^^ = 
1.17 X 10~^^ erg/cm^/s) with fluxes above the detection threshold. 

A similar value, Ncxb ~18, is obtained if we use the source counts in the 0.5-8 keV 
energy band obtained by Kim et al. (2006). We use this number below. 

(2) HMXBs. To estimate the expected number of HMXBs, we use a calibration from 
Grimm et al. (2003) with a modified normalization (Shtykovskiy and Gilfanov 2005): 

A^(> L) = 1.8 X SFR{Lss-^ - 210"°-^), (3) 

where SFR is the star formation rate of the host galaxy. 

According to the IRAS catalog (Rice et al. 1988), the far- infrared luminosity of 
Lfir^I.S X 1O^°L0. The L^r/ii-SFR cahbration (Kennicutt 1998) yields SFRRi 3.1 Mg/yr. 
Using Eq. (3), we then obtain the expected number of HMXBs with luminosities above the 
detection threshold Sii„t(2-10)= 1.17 x lO"-"^^ erg/cm^/s, ^hmxb ~ 19. A slightly larger 
number is obtained if we use the star formation rate derived from the ultraviolet luminosity 
of M51, SFR(UV)Ri 4.3 Mg/yr (Calzetti et al. 2005), which corresponds to ^hmxb ~ 27 
HMXBs. However, this number is less reliable, because the absorption in the ultraviolet is 
uncertain. 

It is worth noting that the calibration by Grimm et al. (2003) also includes the possible 
contribution from SNRs and ultrabright sources, which should be similar to HMXBs in 
spatial distribution. 
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Tablc 1: X-ray scjTirccs tcnvarcl M51 



Source class 


HMXB 


LMXB 


CXB 


Expected number of sources 


19-27 


33 


13-22 


Total observed number 








of sources 




88 





(3) LMXBs. To estimate the number of LMXBs in M51, we use the fact that their number 
is proportional to the stellar mass of the galaxy (Gilfanov 2004). The latter, in turn, is 
proportional to the K-band flux (Bell and de Jong 2001). Using K-band images of M51 from 
the 2MASS Large Galaxy Atlas ( Jarrett et al. 2003) and assuming that the magnitude in a 2' 
aperture is K(< 2') = 6.76 (Rix and Rieke 1993), we obtain K(r< 9.5 kpc)=6.4 for the entire 
galaxy. According to the cahbration by Bell and de Jong (2001) and the color (B-V)o=0.56 
(de Vaucouleurs et al. 1991), this corresponds to M^, = 2.63 x lO^^M©. Using the galactic 
stellar mass - LMXB number calibration (Gilfanov 2004), we then obtain '^lmxb ~33 (22 
in the ring 2.3 kpc<r<9.5 kpc, 13 within the corotation radius 2.3 kpc<r<5.6 kpc) LMXBs 
with luminosities above the threshold value, Sam(2-10)= 1.17 x 10~^^ erg/s/cm^. 

Our estimates roughly agree with the observed total number of sources (sec the table). 
The distribution of X-ray sources in M51 is shown in Fig. 3. It should be kept in mind 
that since the nucleus has a high X-ray brightness, Tcrashima and Wilson (2004) excluded 
the central part of the galaxy from their analysis. Therefore, the central depression in the 
number of sources should not be interpreted as the absence of X-ray binaries in this region. 

The spiral structure in the distribution of X-ray binaries. It is natural to investigate the 
spiral structure in the distribution of X-ray binaries with respect to a certain indicator of the 
region of ongoing star formation or, in other words, an indicator of the shock front location. 
For example, the location of the most massive (and, accordingly, youngest) stars, the inner 
edge of the dust lane clearly seen in the B band, the gas velocity distribution, etc. can 
serve as such an indicator. In our case, it is convenient to use the Hq, emission peak, which 
essentially reflects the location of the youngest stars, as the location of the spiral. For this 
purpose, we used the HST mosaic of M51 (Mutchler et al. 2005). The locations of the spiral 
arms drawn along the emission peak are indicated in Fig. 3 by the solid line. In addition, 
we used bright HII regions (S> 1.5 x 10~^^W/m^) from the catalog by Petit et al. (1996) as 
an additional indicator. 

To estimate the degree of concentration of Xray sources to the spiral arms, we calculated 
the distribution of their distances to the nearest spiral arm using sources at galactocentric 
distances 2.3-9.5 kpc. Positive and negative distances arc assigned to sources above and 
below the spiral, respectively. The derived distribution (corrected for the galaxy inclination) 
is shown in Fig. 4. 

Despite the obvious physical meaning of the constructed distribution, the interpretation 
of its shape is nontrivial. Indeed, a solid angle can be associated with each interval of 
distances on it. As the solid angles for different intervals can differ, the distributions in this 



Fig. 3: HST Hq image of M51 (Mutchler et al. 2005). The circles mark the positions of 
the X-ray sources detected by Chandra (Terashima and Wilson 2004); the circle radius 
reflects the source X-ray luminosity. For a higher contrast, the circles have different 
colors in different regions. Also shown are the two arms of the spiral structure used 
as the location of the region of ongoing star formation. The two circles around the 
galactic center (the light and dark ones at small and large distances, respectively) 
indicate the region from which the sources are taken for our analysis (2.3-9.5 kpc). 
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Fig. 4: (a) Distribution of distances to the nearest spiral arm for X-ray sources (wide 
histogram) and bright HII regions (narrow histogram, reverse hatching). The distribu- 
tion for HII regions is normaUzed to the difference between the total number of X-ray 
sources and the expected number of LMXBs and AGNs. The solid and dotted curves 
correspond to the predictions for LMXBs and AGNs, respectively, (b) The difference 
between the distribution of X-ray sources and the total distribution of LMXBs and 
AGNs (wide histogram) and the distribution of bright HII regions (narrow histogram, 
reverse hatching). 



representation will be distorted, more specifically, they will be more concentrated to zero 
than they actually are. For example, a uniform distribution of sources, such as AGNs, will 
have a maximum near zero whose distinctness will be determined by the shape of the spiral 
and by the boundaries of the region from which the X-ray sources are taken. As a result, 
the observed distribution will have a clear meaning only when it is considered in comparison 
with another, "calibration", distribution. Obviously, in our case, the distribution of bright 
HII regions and a uniform distribution are such distributions. 

The observed distribution is determined by the following factors. 

(1) Finite width of the star-forming region. Even if the star formation in gas clouds 
is assumed to be triggered at a narrow shock front, the star formation process is not an 
instantaneous event, which smears the spatial region of zero-age stars. As was noted above, 
we consider our distribution relative to the distribution of HII regions, which, obviously, 
will reflect the width of the star-forming region. To make a quantitative comparison of the 
distribution for HMXBs with the distribution for HII regions possible, we normalized the 
latter to the difference between the total number of X-ray sources and the expected number 
of LMXBs and AGNs. 

(2) Dynamics of HMXBs. The distortion of the distance distribution due to the displace- 
ment of sources relative to the region of ongoing star formation was demonstrated above. 
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The expected distribution of HMXBs with allowance made for their dynamics is shown in 
Fig. 2. A characteristic feature of the distribution is its asymmetry, which is particularly 
noticeable if only the sources within the corotation radius are taken into account. In ad- 
dition, the distribution of HMXBs is spread due to the initial velocities acquired through 
supernova explosions. 

(3) LMXBs. As was noted above, a high stellar density contrast in M51 {da /a ~2-3, 
mainly within the corotation radius) should lead to a concentration of these sources to the 
spiral arms. Obviously, their spatial distribution should differ from the observed one for 
HMXBs. To estimate the distribution of LMXB distances to the nearest spiral arm, we 
generated a model population of sources distributed over the galaxy in accordance with the 
K-band flux. The number of sources is taken to be large (rs 10^) to suppress the contribution 
from the Poisson noise. The derived distribution for the model sources was normalized in 
accordance with the expected number of LMXBs, ^lmxb — 22. 

(4) Contribution from background sources. To estimate the contribution from background 
sources, we calculate the distribution of distances to the nearest spiral arm for sources 
distributed Tiniformly over the sky and normalize it to the predictions of the source counts 
(see above), N^xs^lS. 

The above components are shown in Fig. 4. Obviously, the distribution of X-ray sources is 
wider and more asymmetric than that of HIT regions. This is the result of a cumulative effect 
from LMXBs, AGNs, and HMXB dynamics. To estimate the contribution from HMXBs, we 
give the difference between the distribution of X-ray sources and the sum of the predicted 
distributions for LMXBs and AGNs. The derived distribution is also asymmetric; it is 
similar in pattern to the prediction of the simplest model for the displacement of HMXBs 
(Fig. 2). However, the result is statistically insignificant. The significance of the fact that 
the distribution of X-ray sources is not a sample from the total distribution of LMXBs, 
AGNs, and bright HIT regions calculated by Kuiper's test is less than 2a (pRi 13%). 

Interestingly, the formal significance of the concentration of X-ray sources to the spiral 
arms is also low, although the reverse seems to be true at first glance at Fig. 3. Indeed, 
judging by Fig. 3, the density of sources near the spiral arms clearly exceeds their density 
in the interarm space. If we do a quantitative count, then we will find that almost as many 
sources are observed within a strip of ±600 pc as outside it, 34 versus 32. At the same time, 
the corresponding areas are ILl versus 19.7 square arcminutes. Thus, the concentration 
is actually present, but its statistical significance is low. The significance of the fact that 
the distribution of distances of X-ray sources to the nearest spiral arm is not a sample from 
the distribution of sources distributed uniformly over the sky calculated by Kuiper's test 
is also found to be low, about 2a for sources with 2.3 kpc<r<9.5 kpc and slightly higher 
for sources within the corotation radius. Thus, the existing number of sources is clearly 
not enough to reach firmconclusions about the distribution of HMXBs relative to the spiral 
arms in M51. Moreover, as we said above, M51 is one of the best examples for studying 
the spiral structure among the galaxies observed by Chandra. Although a comprehensive 
analysis of the archival data on other spiral galaxies is beyond the scope of this paper, we 
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Fig. 5: Approximation of the spiral structure of our Galaxy by a logarithmic spiral 
and a model of th spiral structure from Taylor and Cordes (1993) (thick curve). The 
location of the Solar system and the solid angle containing the tangent to one of the 
spiral arms are shown. 

expect the statistical significance of tlie manifestation of a spiral structure in X-ray sources 
in them to be even lower. Further deep Chandra observations of spiral galaxies with intense 
star formation are required to investigate the manifestation of the spiral structure in the 
distribution of X-ray sources in more detail. 

2.2. High-Mass X-ray Binaries in Our Galaxy 

The simplest kinematic model constructed here also allows us to make predictions 
about the expected distribution of HMXBs in our Galaxy observed from the Solar system 
(r=8.5 kpc). For this purpose, we assume that the spiral structure of our Galaxy follows a 
logarithmic spiral with four arms and a pitch angle if) = 13° that originates at a Galacto- 
centric distance of 3.3 kpc (see, e.g., Vallee 2005, 1995) (Fig. 5). Figure 5 also presents a 
model of the spiral structure in the distribution of free electrons in the Galaxy from Taylor 
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and Cordes (1993). Obviously, its shape deviates from a logarithmic spiral only slightly. We 
assume that the SFR and, accordingly, the number of HMXBs along the spiral arm length 
element decrease with decreasing exponential Galactic disk density. The characteristic width 
of the spiral arm is taken to be w=0.3 kpc (FWHM=0.7 kpc). The SFR at an arbitrary 
point in the Galaxy is calculated as follows: 



where the integral is taken along the spiral arms, dl is the length element of the spiral arm, 
r is its Galactocentric distance, d is the distance from the point to the spiral arm, and 
ro=3.9 kpc (Benjamin et al. 2005). 

After their birth, the HMXBs are displaced relative to the instantaneous locations of 
the spiral arms in accordance with Eq. (2). The rotation curve is assumed to be flat in the 
galactocentric distance range of interest, V/(220 km/s)=ai(i?/8.5 kpc)"^+a3 (ai = 1.00767, 
02 = 0.0394, as = 0.00712) (Brand and Blitz 1993); as the spiral density wave pattern speed, 
we take Qp=24 km/s/kpc (Dias and Lepine 2005), which corresponds to the corotation radius 

~ 9.3 kpc. 

Having integrated the SFR along the line of sights originating from the Solar system, we 
obtained the distribution of HMXBs in Galactic longitude. The distributions for objects 
with ages of 0, 40, and 80 Myr constructed by assuming that we are capable of detecting 
only sources at heliocentric distances smaller than 6.5 and 11.5 kpc are shown in Fig. 6. The 
detection limits roughly correspond to the depths of the INTEGRAL survey for sources with 
luminosities Lx = 10^^ erg/s and 10^^-^ erg/s (Lutovinov et al. 2005). As might be expected, 
the positions of the maxima in the distribution for the youngest objects correspond to the 
spiral arm tangents (see Fig. 5), since in this case the integrated SFR along the line of sight 
is at a maximum. At the same time, as we see from Fig. 5, the maxima in the distribution 
of older objects can be displaced significantly. Besides, additional peaks related to the fact 
that the hitherto invisible inner parts of the spiral arms become visible appear. 

Obviously, the weakness of the constructed model is that we do not know the exact shape 
of the spiral arms and the SFR behavior along them. Indeed, if the SFR is highly nonuniform 
along the spiral arms, then the maxima in the distributions of young sources can also be 
observed in directions different from the spiral arm tangents. 

The displacement of the HMXB maxima relative to the locations of the spiral arm tangents 
were actually observed by Lutovinov et al. (2005) (sec Fig. 7 in this paper, the peak in a 
direction l«i +40°). The observed position of the peak corresponds to binaries with ages 
«i 40 — 80 Myr, which corresponds to the oldest binaries from the standpoint of HMXB 
evolution. However, to reach specific conclusions, we must primarily have a larger number 
of sources and, in addition, understand what contribution the SFR nonuniformity along 
the spiral arms makes to the distribution, hensive analysis of the archival data on other 
spiral galaxies is beyond the scope of this paper, we expect the statistical significance of 
the manifestation of a spiral structure in X-ray sources in them to be even lower. Further 
deep Ghandra observations of spiral galaxies with intense star formation are required to 




(4) 



-14- 




Fig. 6: Model distribution of HMXBs in our Galaxy in longitude for objects with ages 
of (solid curve), 40 (dotted curve), and 80 Myr (dashed-dotted curve). Obviously, 
the solid curve describes the manifestation of the spiral structure in the youngest 
indicators, such as the Hq line and the distribution of the most massive early-type 
stars. We see that the peaks in the distribution of HMXBs are displaced considerably 
as the latter "grow old". The distributions were constructed for sources closer than 
(a) 11.5 and (b) 6.5 kpc. 

investigate the manifestation of the spiral structure in the distribution of X-ray sources in 
more detail. 



CONCLUSIONS 
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We considered the spatial distribution of HMXBs relative to the spiral structure of the 
host galaxy. We constructed the simplest model for the kinematics of HMXBs. It shows 
that they can be displaced to an appreciable distance from the location of the shock front at 
which ongoing star formation takes place over their lifetimes (see Figs. 1 and 2). As a result, 
the spiral structure in the distribution of HMXBs can be displaced relative to the spiral 
structure observed in classical SFR indicators, such as the Hq, emission. The displacement 
should be most pronounced for binaries with neutron stars and minimal for binaries with 
black holes. The evolution of theHMXB number with time elapsed since the star formation 
event can be judged by the displacement pattern. 

As an illustration of this effect, we considered the spatial distribution of X-ray sources in 
the galaxy M51 by analyzing the distribution of their distances to the spiral arms. We showed 
that the components attributable to HMXBs and LMXBs and background AGNs could be 
identified in the distribution using K-band data and background source number counts. The 
distribution of HMXBs shows a clear tendency for them to concentrate to the spiral arms. 
In agreement with predictions of the kinematic model, it is wider and more asymmetric than 
the distribution of bright HII regions, which reflects the region of ongoing star formation 
(see Figs. 2 and 4). However, the statistical significance of both the concentration of X-ray 
sources to the spiral arms and the fact that the distribution of HMXBs differs from that of 
bright HII regions is low (^ 2a). Further deep Chandra observations of spiral galaxies are 
required to reach firm conclusions about the distribution of HMXBs relative to the spiral 
arms. 

We also calculated the expected distribution of HMXBs in our Galaxy observed from 
the Solar system. We showed that the maximum number of sources could be observed in 
directions different from the directions tangential to the spiral arms, where the integrated 

current SFR is at a maximum (Fig. 6), because the HMXBs are displaced relative to the 
shock front. The displacement pattern allows the peculiarities of the HMXB distribution 
observed by Lutovinov et al. (2005) to be explained qualitatively. 
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